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On the basis of the Ward-Pitaevskii identity, the residual resistivity po is shown to exhibit 
huge enhancement around the quantum critical point of valence transition in Ce-based heavy 
electron systems. This explains a sharp peak of po observed in CeCu2Ge2 under the pressure at 
P --^IGGPa where the superconducting trasition temperature also exhibit the sharp peak. 



I 

o 
o 



> 

00 

in 



o 



o 
o 



X 
J3 



KEYWORDS: enhanced impurity scattering, critical valence fluctuations, CeCu2Ge2, CeCu2Si2 



Recently, renormalization effect of impurity potential 
by quantum critical fluctuations has begun to attract 
much, attention, while the effect of impurities on 
the universality class of critical fluctuations was clari- 
fled quite long ago,*^ and that on the temperature de- 
pendence of the resistivity at quantum criticality has 
been discussed recently.^' A theoretical guideline for dis- 
cussing exactly such an effect has already been put forth 
more than three decades ago by Betbder-Matibet and 
Nozieres in the framework of the Fermi liquid theory. ^"^^ 
Indeed, they showed on the basis of the Ward identy that 
the impurity potenital, in one-component Fermi liquid, 
is renormalized by many-body effect in the forward scat- 
tering limit as 



u{k ^ 0) 



1 



z(l + 



-u{k 0), 



(1) 



where u(fc) is the bare non-magnetic impurity potential 
and u{k) are the renormalized one, z is the renormaliza- 
tion amplitude including all the manybody effects, and 
Fq the Landau parameter relevant to the correction of 
the charge susceptibility. For the potential of magnetic 
impurity, the relation similar to (1) holds with the Fermi 
liquid parameter Fq being replaced by Fq which gives 
the Fermi liquid correction of the spin susceptibility. 

Heavy electron compound CeCu2Ge2 at ambient pres- 
sure changes drastically its electronic state at the pres- 
sure P ~ 17GPa where the coefflcient A of T^-term of the 
resistivity p{T) decreases by about three orders of mag- 
nitude and the universal ratio 7 being the Som- 
merfeld constant, changes from the value of heavy elec- 
torns to that of conventional c?-band metals decreasing by 
25 times. This suggests that the rapid valence change 
of Ce ion occurs at around that pressure. Indeed, the 
rapid volume change maintaining the crystal symmetry 
was observed at around the "critical" pressure by SOR- 
X ray diffraction implying that the rapid valence change 
occurs there. It was also observed that the residual 
resistivity po exhibits sharp peak at the same pressure. 
Similar behavior has been observed in CeRhlns recently 
discovered pressure induce superconductor.^^) Recently, 



it was shown that po can be enhanced through the renor- 
malization of impurity potential by exchanging the crit- 
ical valence fluctuations."^'"'' However, the argument 
is based on the perturbational treatment with the use of 
the phenomenological form of valence fluctuation prop- 
agator. A purpose of this Letter is to extend that 
treatment so as to take into account the full effect of 
vertex corrections to the impurity potenital on the basis 
of the Ward-Pitaevskii identity in two-component Fermi 
liquid. Namely, we derive a technically exact formula for 
the renormaliztion of impurity potential for the forward 
scattering associated with a rapid critical valence change 
in heavy electrons form the Kondo regime to the valence 
fluctuation one. 

We start with an extended periodic Anderson model 
(PAM), 



H 



per 



P^paJpa 



(2) 



where the conventional notations for PAM are used ex- 
cept for C/fc, the f-c Coulomb repulsion. It has recently 
been shown that the effect of J7fc is important for the 
rapid valence change to occur as the f-level Cf is in- 
creased approaching the Fermi level by the effect of pres- 
sure.^'*' The one-particle Green function for a given 
spin direction in this system is given formally as 

Gs{p,e) Gfc{p,e) ~^ 
Gct(p,e) Gcc(p,e) 



[G-.i(p,£)] = 



(3) 



-V* 



-p - S;ff(p, e) 



-Vp- Sfc(p,e) 
£ - £.P - ^ccip, e) 



(4) 



where Eg is the selfenergy of f-electron and there also 
exist Efc and Ecc as the many-body effect due to Ug 
and C/fc. It is noted that G^^^ means zj-element of the 
inverse matrix of the Green function, e.g. Gii =/= Gg^ 
while Gil = Gff. 

For coupling to nonmagnetic impurities, the scatter- 
ing matrix for this system is given generally as Uij{k) = 
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Sa=i '^a,{k)Xfj where ui and U2 are the variations of po- 
tential on f-electrons and c-electrons, respectively, while 
U3 and M4 represent the strength of the f-c mixing scat- 
tering. Here Xfj is the bare vertex for coupling to impuri- 
ties and is given by [A^ = (l + r^)/2, [A^ = (1-t^)/2, 
[A3-] = r^/y2 and [Xjj] =Ty/V2 with f being the Pauh 
matrices in the f-c space. 

In the following, a Greek index, e.g. a, represents 
the dependence on both the component i = 1,2 and the 
spin (T =t, I, and the summation is assumed to be taken 
for repeated indices. Then, the one-particle Green func- 
tion and the scattering matrix mean the tensor product 
multiplied by the unit matrix with respect to the spin 
variables. 

The Ward-Pitaevskii identity relevant to the present 
problem is given by considering the linear response of 
Grya caused by the shift of the parameter denoted by ea 
with the chemical potential jj, being fixed. Here Sa is the 
f-level et, the center of the conduction band Ec or f-c mix- 
ing V, i.e., *[ea] = (£1,62,63,64) = (6f , 6c, ^2^' , ^21/") 
with V and V" being the real and imaginary part of V, 
respectively. One can show, by anlyzing the structure of 
perturbation series of the selfenergy, that the following 
identity holds: 



dea 



\ dea 



(5) 



(2^)- 



where F*^ is the so-called /c-limit of the full vertex func- 
tion, and {GG}k is the same limit of particle-hole Green 
function pair with the four-vector abbreviations p — 
{p,s), etc. The process of renormalization of impu- 
rity potential is represented by the Feynman diagram 
as shown in Fig. 1. In the limit of forward scattering, 
i.e., fc 0, the renormalized potential Ua{k) and the 
bare one Ua{k) are in the relation 

lim Ua{k) = lim ^A" -U6(fc) 



fe^O 



2 



(2^)' 



.(7) 



The reason why the fc-limit vertex appears in (7) is that 
the impurity scattering is elastic maintaining the energy 
transfer w = 0. Therefore, by the relation (6), the impu- 
rity potential is renormalized by —{dG~^{p,e)/dea)fi in 
the limit k ^ 0. 

In heavy electrons, the important effect of impurity on 
the quasiparticles arises from the variations of potential 
on f-electrons, because the quasiparticles consist mainly 
of f-electrons. Of these effects, those by displacement 
of non-f elements from the regular alignment around Ce 
ions is subject to remarkable renormalization by the crit- 
ical valence fluctuations, because the impurity potential 
due to such effects is off the unitarity limit and has a 
space to be renormalized furthermore. On the other 
hand, the defect of Ce ions gives rise to the unitarity 
scattering from the beginning in heavy electron state so 
that its potential is subject only to a gradual renormal- 



ization with a weak anomaly around a possible transition 
point from the Kondo regime to VF one. It is noted that 
even if we consider only the effect of the shift of the 
f-level, i.e. in the case of U2 — — = 0, the many- 
body effect due to Ufc gives rise to the effective f-c and 
c-c scattering as can be seen in (7). 



lim 



iiAk) 





Fig. 1. Foynman diagram for the exact vertex correction of im- 
purity potential u(k). 



In the present system described by the Hamiltonian 
(2), the total number density of f-electrons and conduc- 
tion electrons is conserved. This leads to the following 
identity: 

dG-l{p) _ ^ 



de 



d^g 

(2^ 



dY.-fa{p) 

de 



(8) 



r75,./3(p,9){G^c(9)Gc5(g)}. (9) 



where F" is the so-called w-limit of the full vertex func- 
tion and {GG}i^ is the same limit of particle-hole Green 
function pair. The relation between F*^ and F" is given 
by 

X [{Ga{(l)G^.{q)}k - {G^^{q)G^MUT\,^^p{q,p')- (10) 

By using formulas (9) and (10), we can evaluate the left- 
hand side of (5) as shown below. 

First, we write the total number density n in terms of 
the Green function as 

d^p 



{2^y 



-Gaa{p)- 



(11) 



Differentiation of (11) with respect to with ^ being 
fixed gives 

r^{Gaf}{p)Gp^{p)}k \ - 



dn 
dea 



{2^Y 



dea 



(12) 



By (6) and (10), we find the right hand side of (12) is 



dV 
(2^ 



d^q 



i2n) 



j{Gat3{q)G,}y{q)}^T':^^^^^{q,p) 



X [{Ga{p)G,Ap)}k - {Gcdp)G,.{p)}^ 



fdG^ 
[ dea 



d^p 



SkC - 



d^q 
(2^ 



(27r)4 

x{Ga{p)Gr,4p)}^^1, 

According to the Ward identity (9), the integrand of the 



{Ga(3 {q)G(3^ (q) }a;F^^_„ J (g, p) 

(13) 
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second term in (13) is just X^rjdGr]^{p, e)/de so that this 
term equals zero. Substituting (9) into the first term in 
(13), we obtain 



dn 
dea 



V dea 



X [{G0^ip)Gsc.{p)}k - {Gfs^ip)Gsc.ip)U . (14) 



Near the Fermi level, (3) is expressed in terms of quasi- 
particle as 



G7^\p, s) « a7^\p)e - eijip), (15) 
= ^Hp) = 9Gr.'^{p,e)/de\s=o and = 



where a- = 
^ij{p ) is given by 



■ (16) 



ef-/x + I]ff(p,0) y + Sfe(p,0) 
y*+Srf(p,0) ^jr+Eec(p,0) 

Then the one-particle Green function G can be written 
in the form 



Gij {p, e) = Aij {p)/[e- E^] + nonsingular part, (17) 

where the renormalization amplitude Aij (p) and the dis- 
persion of quasiparticles Ep near the Fremi level is given 

by 



[Aij{p )] = det(a) /tr(ae) 



^p- + Scc(p,o) -y-Sfe(p,o) 

-F*-Sef(p,0) ef-M + Sff(p,0) 



Ei^- 

P 2 



tr (ae) ± ^ [tr (ae)] 
respectively, with a = [a^ ] and e 



4det (ae) 



,(18) 



(19) 



[dj]. By(17),thedif- 



fernce between the fc-limit and the w-limit of the product 
GG can be obtained to be 

{Gij{p)Glm{p)}k - {Gij{p)Glm{p)}oj 

= -2wiAij{kF)Aim{kF)S{e)6{EZ). (20) 

Substituting this equation into (14) and taking the sum- 
mation of the spin variables, we obtain 



fdn\ ~ ~ /aG-/(fcF,o) 
(^j^ = iVp^.(M 



(21) 



where iVp is the density of states of the quasi-particles. 
In deriving (21), we have used the relation det[>ly(A;F)] — 
so that Aij{kF)a~j^ {kF)Aim{kF) = Aimikp). By virtue 
of (7), wc finally find that the renormalized potential 
acting on the quasi-particles is given by 

dn 



A\kF)ubik 



Nf \dea 



Ua{k ^ 0), (22) 



where A''{kF) = X'ljAji{kF)- 

Equations (21) or (22) can be obtaind more easily 
in the limiting case of heavy electrons, af = an <C 1, 
ai2 = 0,21 ~ and 022 ~ 1, which is the only way for 
the quasiparticlc to acquire the extremely heavy mass. 
In such a case, the quasiparticles consist mainly of f- 
electrons whose weight in the quasiparticlc state is given 
by 1 — ctf|l^P/Cfep ~ 1- Therefore, the derivative with 



respect to ef in (5) is approximated in a high accuracy, 
neglecting af compared to 1, as 



1 s 



£=0 



/ dkF 



(23) 



where the velocity of quasiparticles is defined as Vp = 
dE~/dp, and p and e have been approximated by p = 
fcp and £ = 0, respectively. In deriving (23), wc have 
considered that the dispersion of quasiparticles near the 
Fermi level is approximated as E^ !v v~{p — kp) and 



p— >/cF \ 0€f 



fdkp 



■ (24) 



Since the renormalization factor Of included in and 
that in the denominator of (23) cancels with each other, 
we obtain 



d_ 
dcf 



Gi\p,e) 



1 



dkp 
dnf\ 



(25) 
(26) 



where Uf is the f-electron number density and Nf is the 
density of states of non-interacting electrons described by 
(2) at the Fermi level. The reason (25) is approximated 
by (26) is as follows: A variation of fcp under n being 
fixed corresponds mainly to that of f-electron number, 
because the quasiparticles consist mainly of f-electrons-^'') 
and the change of conduction-electron number is limited 
by the fixed chemical potential /x. This physical picture 
can be verified on the basis of Gutzwiller approximation 
appUcd to PAM without Ufc-^^'^^^ 

By the relations (6), (7), and (26), the impurity po- 
tential acting on f-electrons, in the forward scattering 
limit, is renormalized by the valence fluctuations associ- 
ated with the crossover from Kondo regime to VF one 
as 



u{k 0) 



1 / dnf 



u{k 0). 



(27) 



This is an anlog of the relation, (1), in which the factor 
l/z(l + Fq) is reexpressed as Xchargo/-^F- The relation 
(27) implies that the impurity scattering is critically en- 
hanced if the valence of Ce-ion changes critically as the 
f-level Cf is tuned, relative to the Fermi level, by the pres- 
sures. Indeed, it has been demonstrated theoretically 
that the derivative — (9nf/9ef)„ can diverge in the sys- 
tem described by the model Hamiltonian (2) with appro- 
priate values of Ufc and ef This means — {dnf / def)^ 
also diverges there, because the following relation holds, 
up to the approximation (26), 



dnf 
dcf 



duf 

d^ 



dnf 
dn 



(28) 



where the derivative —{dnf/dn)^^ is a small number of 
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the order of z, the renormaUzation ampUtude. 

In order to see how this enhancement of impurity po- 
tential affects the behaviors of the resistivity, we need to 
know the fc-dependence of u{k) for the scattering from 
p— k/2 to p + fc/2 near the Fermi surface. Although it 
is not easy to determine the ^-dependence accurately in 
general, it may be reasonbale to parameterize as 

^^^^ " ;^tW"(^)' 

where r] is the inverse valence susceptibility, ry"^ = 
|(5nf/5ef)^|/7VF, and ~ In the case where 

the bare impurity potential causes essentially the Born 
scattering, the enhancement of the residual resistivity po 
by the critical fluctuations becomes gigantic. Indeed, po 
is given as 



27r7VFCimp|w(fc)|^(l - COS 

[r? + A(2fcFsin^(6l/2)] 




(30) 



where Cimp is a concentration of impurity, 6 is an angle 
between p±k/2, and the on-shell condition £=^-^^^2=0 
has been used. (• • •)fs nieans that the average with re- 
spect to p over the Fermi surface is taken. Here it is 
noted that explicit dependence of renormaUzation ampli- 
tude does not appear due to cancellation between that for 
DOS and that for the damping rate of quasiparticles.^"-* 
Calculation of angular average over 6 is performed easily 
giving rise to 



po (X In 



1 

v' 



(31) 



This result remains valid even if we take into account 
higher order terms by calculating the t-matrix. It is be- 
cause the scattering by the renormalized potential (27) 
is nothing but the Rutherford scattering in the limit of 

It should be noted that the present result is not con- 
tradictory to the Friedel sum rule according to which the 
scattering probability does not diverge so long as the ex- 
tra charge accumulated locally around the impurity is 
finite. The form of renormalized impurity potential 
(27) becomes long ranged as ^ even though the bare 
potential is short ranged. Namely, the change of valence 
near the impurity site extends in long range proportion- 
ally to 1/r, r being the distance from the impurity site. 
As a result, total amount of valence change around the 
impurity from that of the host metal is divergent while 
the local charge of f- and conduction electrons remains 
finite. So, the effect of impurity becomes long ranged 
making the scattering probability divergent. 

The expression of po given by (31) can explain a huge 
enhancement observed in CeCu2Ge2 and CeCu2Si2 at 
around the critical pressure where the rapid valence 
change seems to occur. This huge enhancement should 
be compared to the moderate enhancement around the 
magnetic quantum critical point where the enhancement 
arises only through the renormaUzation amplitude z as 
discussed in Ref.''' TUc present mecUanism of huge en- 
hancement of pq is related to other systems near tUc 
quantum critical point associated with charge instabil- 



ity. For example, such a huge enhancement has been 

observed in Cd2Re2 07 at around the pressure where 
the charge ordering temperature vanishies and the su- 
perconducting transition temperature Tc is considerably 
enhanced compared to that at the ambient pressure. 
Similar enhancement is expected in /3NaVgOi5 which ex- 
hibits pressure induced superconductivity of T^. — 10 at 
around the pressure where the charge ordering is sup- 
pressed. 

The enhancement of the two-dimensional charge suc- 
septibility was observed by pUotocmission spectroscopy 
near tUe metal-insulator transition in Uigli-Tc ciiprates 
LSCO.^'"'^ This opens a new possibility tUat the doped 
ions, Sr ions in LSCO, which have only subtle influence 
on electrons in Cu02 plane in the over and optimum re- 
gions, are transformed into the strong scattering center 
in the under doped region. Therefore, the carrier doping 
can suppress of d-wave superconductivity consider- 
ably near the metal- insulator phase boundary. 
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